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ABSTRACT 



X 



We have analyzed the optical and radio properties of a sample of 3CR FR I radio galax- 
ies which have Hubble Space Telescope imaging capable of detecting optical cores. The jet 
powers of the FR I radio galaxies are estimated from their low-frequency radio luminosities, 
and the optical core luminosity is taken as an upper limit on the emission from any unob- 
scured accretion disc. We argue that if the accretion discs in these sources are assumed to be 
advection dominated accretion flows (ADAFs), or adiabatic inflow-outflow solution (ADIOS) 
flows, then the Blandford-Znajek mechanism ( iBlandford & Znaiei3ll977^ provides insuffi- 
cient power to explain the high radio luminosities of at least a third, and perhaps all, of the 
sample. We suggest instead that a significant fraction (the 'high-jet-power' third), and per- 
haps most, of the 3CR FR I radio galaxies have normal accretion discs, but that their optical 
cores can be hidden, with any HST-detected optical synchrotron emission coming from jets 
on scales larger than the obscuring material. A normal accretion disc hypothesis, at least for 
the high-jet-power third of the 3CR FR Is, explains why narrow-line luminosity correlates 
with radio luminosity. It also explains why one object in the sample (3C 386) has an observed 
broad-line nucleus. We conclude that there is no evidence to suggest that there is a difference 
in accretion mode between FR I and FR II radio galaxies. 

Key words: accretion, accretion discs-galaxies: jets-galaxies: nuclei 



• ^ 1 INTRODUCTION 



FR I radio galaxies (defined by edge-darkened radio structure) 
have lower radio power than FR II galaxies (defined by edge- 
brightened radio structur e due to compact jet-terminating hot spots) 
jFanaroff&Rile\ll974 . What causes the difference between FR I 
and FR II radio galaxies is still unclear. The different explanations 
of the division of FR I and FR II radio galaxies invoke either the in- 
teraction of the jet with the ambient medium or the intrinsic nuclei 
properties of accretion and jet formation processes (e.g., Bicknell 
1995; Gopal-Kiishna & Wiita 1988; Reynolds et al. 1996; Baum, 
Zirbel & O'Dea 1995; Blundell & Rawlings 2001). In the frame 
of unification schemes of active galactic nuclei (AGNs), FR I ra- 
dio galaxies are believed to be the misaligned BL Lac objects, and 
FR II radio galaxies correspond to misaligned flat-spectrum radio 
quasars (see Urry & Padovani 1995 for a review). The central struc- 
ture and activity of radio galaxies are important for understanding 
the physics at work in their nuclei. 

Recently, the optical nuclei of low-redshift radio galax- 
ies have been observed by the Hubble Space Telescope 
(HST), and central compact optical cores detected in many 
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cases ^Caoetti & Celotti"l999';'Chiaberse. Caoetti. & Celotti'l99g|; 
[ Chiaber ge. Capetti & Celotti 2002; Verdoes Kleiin et al. 200^. 
The fact that the optical fluxes of the central compact cores in 3CR 
FR I radio galaxies correlate linearly with their radio core fluxes, 
and the lack of scatter in this correlation, argues for a synchrotron 
origin for the optical emis sion of the cores in FR I radio galax- 
ies i Chiaberge. Capetti. & Celotti 1999) . FR II galaxies, typically, 
show brighter optical cores than FR I galaxies and a more com- 
plex correlation between radio and optical core luminosities sug- 
gests that the disc emission may cont ribute a significant fractio n 
of th eir optical core fluxes (Chiaberge. Capetti & Celottill2002h . 
jChiaberge. Capetti. & Celotti il999.) estimated the upper limits on 
the disc emission of FR I galaxies from the observed optical core 
fluxes and found that the optical core luminosities of FR I radio 
galaxies are far lower than the Eddington case (I/Edd), if the black 
holes in FR I radio galaxies have masses around lO' Mq ; typically 
AI/c,opt/iEdd < 10"", AI/c,opt /iEdd — 10"*" in the case of the 
archetypal nearby FR I galaxy M87. 

In the frame of unification schemes, the viewing angle of 
FR I galaxies is, on average, large, and it is therefore to be ex- 
pected that any geometrically and optically thick torus present 
could hide the nucleus and broad-line emission in most cases 
jUrrv&Padovanilll995h . In the receding torus model iLawrenc^ 
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Il99ll) . in particular, tiie optical emission from a comparatively 
weak disc is likely to be hidden at all but very small angles 
to the line of sight, as the t orus is predicted to lie extremely 
close to the nucleus jffill. Go odrich & DePoJ Il99d) . However, 
IChiaberee. Capetti. & Celottil jl999) argued that high detection 
rate of optical central compact cores of FR I radio galaxies (^85%) 
in the 3CR catalogue suggests such a torus is not present in FR I 
galaxies. For example, the linear resolution achieved by their HST 
observations at the distance of M87 corresponds to ^^6 pc, implying 
that, in this source at least, the synchrotron radio and optical core 
emission is not being hidden by a torus, so if one is present it must 
be less than 6 pc in size. Constraints on the size of any obscuring 
torus clearly worsen, however, for the more distant, powerful FR I 
radio galaxies in the 3CR sample. 

The division between FR I and FR II radio galaxies is clearly 
shown by a line in the plane of the opti cal lumi nosity of t he hos t 
galaxy and the total radio luminosity bv lLedlow & OwenI il996l) . 
The FR I/II division line is a function of the optical luminosity of 
the host galaxy (roughly proportional to the optical luminosity). 
The total radio luminosity of almost all FR II radio galaxies lies 
above the division line, while FR I radio galaxies are in the region 
below the line. They suggested that the accretion rate may play an 
important role in the division between FR I and FR II galaxies. 
Recently, Ghisellini & Celotti (2001) used the optical luminosity 
of the host galaxy to estimate the mass of its central black hole, and 
used the radio luminosity to derive the jet power. They proposed 
that the FR I-FR II separation can be interpreted as occurring at 
a constant ratio between the jet power and the black hole mass. 
This suggests that the FRI/FRII division is linked to physics on 
very small scales, and potentially suggests a difference in accretion 
mode. 

For radio galaxies and quasars, their optical/UV continuum is 
a mixture of the emission from the jet and from the disc, which 
prevents estimating the emission from the disc directly from opti- 
calAJV continuum. Moreover, obscuration of the accretion disc is 
known to happen in some radio galaxies jSimpson. Ward & Walj 
I2OOOI) so an orientation-independent estimator of the nuclear lumi- 
nosity is needed. 

A widely adopted approach to estimate the accretion disc 
emission in radio galaxies and quasars is to use their narrow emis- 
sion lines (e.g., Rawlings & Saunders 1991; Xu, Livio & Baum 
1999), or broad lines (e.g., Celotti, Padovani & Ghisellini 1997; 
Cao & Jiang 1999; 2001), which are ionized by the photons emit- 
ted from the disc (although photo-ionization by radiative shocks 
driven by the radio source can also play an important role: e.g., 
Inskip et al. 2002). The applicability of this approach to FR I ra- 
dio galaxies has been questioned (Zirbel & Baum 1995) but it is 
at least roughly in accord with the correlations between line and 
optical luminosities seen in high-spatial-resolution observations of 
nearby FR I radio galaxies (Verdoes Kleijn et al. 2002). 

The cosmological parameters Ho = 50 kms~^ Mpc~^ and 
qo = 0.5 have been adopted in this paper. This is not now the con- 
ventional choice, but as we are studying only low-redshift objects, 
this has no significant bearing on the results of this paper. 



2 SAMPLE 

IChiaberge. Capetti. & Celottil < 1 999^ processed the data of HST ob- 
servations in the public archive for 32 of all 33 FR I radio galaxies 
belonging to the 3CR catalogue (3C76.I has not been observed). 
In this work, we use these 33 FR I radio galaxies belonging to 



the 3CR catalo gue as our sample. Key data on the se sources are 
listed in Table I.'Chiaberge. Capetti & Celotti" (■2002^ detected cen- 
tral compact cores unresolved at the HST resolution in 22 radio 
FR I galaxies. Five sources show larger optical cores (> 0.15"), 
which are different from those sources with detected central com- 
pact cores (0.05" — 0.08"). For these five sources, only upper limits 
on the central compact core emission are given. The remaining five 
sources have complex nuclear optical structures. In four of these 
five sources, bright compact knots are resolved and offset from the 
center of the galaxy. The optical image of 3C315 has an elongated 
structure and no central source is seen. No reliable optical core flux 
is available for these five sources. 



3 CENTRAL BLACK HOLE MASS 

In order to estimate the central black hole masses of FR I galaxies, 
we use the relation between host galaxy absolut e magnitude Mr at 
J?-ban d and black hole mass Mbh proposed bv lMcLure & Dunlopl 
J2002I) 

logio(Mbh/M0) = -0.50(±0.02)Mfl - 2.96(±0.48). (1) 

The derived black hole masses of all sources in the sample are listed 
in Table 1. For M87(3C274), i ts central black hole mass has been 
estimated to be 2:^ 3 x lO^Mp, jFord et alJl99llHarms et alJl994 
lMacchettoetaIlI997^ . 



4 JET POWER AND ITS ORIGIN 

The relation betwe en jet power and radio luminosity proposed by 

Ojet ^ 3 X 10^V'/'^?^I W, (2) 

where L151 is the t otal radio luminosit y at 151 MHz in units of 
10^** W Hz"^ sr~ MWiIlottetalJ<I999 () have argued that the nor- 
malization is very uncertain and introduced the factor / to account 
for these uncertainties. T hey use a wide v ariety of arguments to 
suggest that 1 < / < 20. lBlundell & Rawlings I.20QQ.) argued that 
/ 10 is a likely consequence of the evolution of magnetic field 
strengths as radio sources evolve. This would mean that the bulk 
kinetic (jet) and radiative (quasar) outputs of radio-loud AGN are 
of similar magnitude (Fig. 7 of Willott et al. 1999). This relation 
was proposed for radio FR II galaxies and quasars. We adopt this 
relation to estimate the power of jets in FR I radio galaxies which 
is an approximation (because the efficiency factor / for FR Is may 
well be systematically different), but one that is probably good to 
an order of magnitude. For example, we use this relation to esti- 
mate the power of the jet in M87, we find Qjot ~ 6 x 10^^ W 
in the case of / = 1, wh ich, for reasonable / is con sistent 
with previous estimates (e.g., Bicknell & Begelman found 
-,^q36-37 ^. Young, Wilson & Mundell |2002) derived the temper- 
ature and pressure of the gas in its inner radio lobes (^ 1.5 kpc) 
from the Chandra X-ray observation. The cavity of the X-ray gas 
is assumed to be created by the jet and its age is estimated from 
the free-fall timescale. The minimum power of the jet in M87 is 
therefore estimated as ~ 3 x 10^^ W). 

It is evident that the accretion rates of FR I galaxies are very 
low. It has been suggested, therefore, that their accretion discs are 
probably in an advection d ominated accretion flow (ADAF) state 
jGhisellini & Ce lotti"2001). An ADAF radiates inefficiently and is 
very hot. Most gravitational energy of the accreting matter released 
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is carried into tlie black liole. Unlike the optically thick accretion 
disc case, the re is a theoretical upper limit on the accretion rate 
for an ADAF JNaravan & Y i 1995). The optically thin ADAF will 
not exist and it will transit to an optically thick disc, if its accre- 
tion rate is greater than the critical one rhcrit = Merit /M-^dd (see 
Narayan 2002 for a review). The exact value of the critical accre- 
tion rate rhcrit is still unclear, depending on the value of the disc 
viscosity parameter a, i.e ., merit — . 28q^ iMahadevan, . 1 997i) . 
iGhisellini & Celottil i200ll) estimated the power of the jets in ra- 
dio galaxies from their radio luminosity. They proposed that the 
division between FR I and FR II radio galaxies can be expressed as 
a constant accretion rate between ~ 10"^-10~'^ of the Eddington 
one on the assumption of a constant conversion of the jet power to 
bolometric luminosity. This is consistent with the critical accretion 
rate expected by ADAF models in the sense that FR I radio galax- 
ies could be radio sources whose central engines are in the ADAF 
state, and FR lis could be radio sources whose central engines are 
standard, optically- thick accretion discs. It is already known that 
there would have to be exceptions to any such rule: 3C386, an FR I 
in our sample (Table 1), has a relatively bright compact optical core 
and a broad-line nucleus (Simpson et al. 1996); there is at least one 
example of an optically bright quasar with FRI-like radio structure 
(Blundell & Rawlings 2001). 

If the black holes in FR I galaxies are spinning rapidly, the 
rotational energy of the black holes is expected to be transferred 
to the jets by the magnetic fields threading the holes, na mely, 
the Blandford-Znajek mechanism ( Blandford & Zna iekll977h . The 
power extracted from a maximal spinning black hole was calcu- 
lated by[Armitaee & Nataraian ( 1999). Here, we calculate the jet 
power in a similar way. The power extracted from a spinning black 
hole is given by (e.g. Ghosh & Abramowicz 1997; Macdonald & 
Thorne 1982) 

Lbz = ■^ujpB'iRlca^, (3) 

for a black hole of mass Afbh and dimensionless angular momen- 
tum a, with a magnetic field B± normal to the horizon at i?h. _Rh 
is proportional to the black hole mass Mbh and is a function of a. 
Here the factor cop = Qf{Qi^ — SIf)/^^ depends on the angular 
velocity of field lines Sip relative to that of the hole, flh. In order to 
estimate the maximal power extracted from a spinning black hole, 
we adopt cup = 1 /2. 

The magnetic field threading the black hole is maintained 
by the currents in the accretion disc surrounding the hole 
jOhosh & Abramowic z 1997; Livio. Ogilvie & Pringle 199^. The 
pressure of an ADAF is given by ( .Narayan & YL1995.) 



2(5 + 2e') , „ 
C3 = — g{a, e ^ 



p = 1.71 X lO'^a-^c- 



(4) 



where the dimensionless quantities are defined as follows: 

_ Mbh _ Rc^ . _ M 
"""M^' '"-2GM^,' """m^' 



and 



Me 



1.39 X 10^'m kgs"\ 



(5) 



Here the standard value of the accretion ef ficiency factor rj^g = 0.1 
is adopted. The parameters ci and C3 are jNaravan & Yil 19951) 



(5 + 26') , 
ci = — — — g{a, e ) 



3a2 



(6) 



9a2 



where 



f 5/3 - 



/adv V 7 - 1 



and 



1 1/2 



(5 + 2e')2 



1. 



(7) 



(8) 



(9) 



and 



The parameter /adv, which lies in the range 0—1, is the fraction of 
viscously dissipated energy which is advected; 7 is the ratio of spe- 
cific heats. So, the value of the parameter e' is in the range of — 1. 
As done by Armitase & Nataraian ( 1999), we assume B± ~ B, 
Pmag ~ -B^/Stt ~ ap, and two parameters a — 1 and € — 1 
are adopted to maximize the pressure (see Eq.|4j. It is obvious that 
the magnetic pressure is proportional to the dimensionless accre- 
tion rate m (see Eq.|4}. As the accretion rate of an ADAF should 
be less than the critical one merit, we can calculate the maximal 
jet power extracted from a spinning black hole of mass Afbii and 
dimensionless angular momentum a from Eq. J3}' if m = rhcrit is 
substituted into Eq. J4j. T he results of the extreme case for a = 1 
have already been given in lArmitage & Nat araian ('l999|). 

The ADAF solutions may be modified to adiabatic inflow- 
outflow solutions (ADIOSs), if a powerful wind is present to carry 
away mass, angular mo mentum and energy from the accreting gas 
(Blandford & Begelmai] fT999l) . In this case, the accretion rate of 
the disc is a function of radius r instead of a constant accretion 
rate along r for a pure ADAF. For ADIOSs, the gas swallowed 
by the black hole is only a small fraction of the rate at which it 
is supplied, as most of the gas is carried away in the wind before 
it reaches the black hole. The accretion rates at the inner edge of 
the disc for ADIOSs would be at least as low as that required by 
pure ADAF solu tions, and they have sim ilar structure at the inner 
edge of the disc ilChang. Choi & Yil2002h . The jet power extracted 
from a spinning black hole by the Blandford-Znajek mechanism for 
ADIOS cases is therefore similar to that for pure ADAF cases, if 
they have similar accretion rates at the inner edge of the disc. 



5 RESULTS AND DISCUSSION 

The derived jet powers and black hole masses for 33 FR I galaxies 
are listed in Table 1. We searched the literature for flux data on 
the narrow emission lines Hci-l-[Nll]. For those where the flux of 
Ha+[N II] is unavailable, we use [O III] flux multiplying a factor 
0.83 to convert to Ha+[N II] flux as done by Zirbel & Baum ( 199^) . 
We found 25 sources with narrow-line emission data. 



5.1 Jet power versus black hole mass 

As discussed in Sect. 4, the uncertainties in the jet power estimate 
are included in the factor / in Eq. J2j. We plot the relation between 
the jet power Qjct and black hole A/bh in Fig.Qfor / = 1. As dis- 
cussed in Sect. 4, the maximal jet power extracted from a spinning 
black hole with Afbii and a can be calculated for an ADAF sur- 
rounding the hole. The three-dimensional M HD simulations sug- 
gest the viscosi ty a in the discs to be ~ 0. 1 jArmitagd[l993) . or 
- 0.05 - 0.2 taawlev & Balbusll2002h . We adopt m^rit = 0.01 
in our calculations, corresponding to a ~ 0.2, which is a conser- 
vative choice. The calculated maximal jet power is plotted in Fig. 
Q We find that the jet power for more than one-third of the sources 
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Table 1. Data of the sources. 



Source 


z 






' ' ''±1 


Refs. 


1 Off 1/1 ^J^■ 
olO line 


Refs. 


login Oi^t//'^/^ 


login A/kk /Mri 
olO bh/ 


(V) 


(2) 


(■31 


(A) 




(61 


(7) 


(81 


(91 


(101 


3C 28^ 


0.1953 


<30.76 


<22.48 


18.34 


d96 


35.07 


ZB95 


37.1 


8.1 


3C 29 


0.045 


30.71 


23.90 


13.02 


GOO 


33.69 


ZB95 


36.0 


9.1 


3C31 


0.0169 


30.27 


23.05 


11.92 


S73 


33.18 


V02 


35.3 


8.6 


3C 66B 


0.0213 


30.99 


23.55 


12.55 


M99 


33.64 


V99 


35.8 


8.5 


3C 75 


0.0232 




22.95 


11.88 


S73 


32.64 


ZB95 


35.7 


9.0 


3C 76.1 


0.0325 


not observed 


22.55 


14.06 


S73 






35.6 


8.2 


3C78 


0.0287 


31.93 


24.53 


12.50 


M99 


34.02 


ZB95 


35.7 


8.9 


3C83.1 


0.0251 


29.59 


22.75 


12.39 


S73 






35.8 


8.8 


3C 84 


0.0176 


32.31 


25.75 


11.03 


G77 


35.21 


ZB95 


35.9 


9.1 


3C 89=^ 


0.1386 


<30.25 


24.58 


16.55 


d96 


34.68 


ZB95 


37.0 


8.6 


3C 264 


0.0217 


31.37 


23.60 


13.03 


M99 


33.39 


ZB95 


35.6 


8.3 


3C 270 


0.0075 


29.09 


22.87 


10.16 


S73 


<32.46 


ZB95 


35.1 


8.6 


3C 272.1 


0.00334 


29.45 


21.94 


9.26 


S73 


32.67 


ZB95 


34.1 


8.2 


3C 274 


0.0044 


30.51 


23.51 


9.03 


S73 


33.77 


ZB95 


35.8 


9.5 


3C 277.3=^ 


0.0853 


30.69 


23.58 


15.71 


E85 


34.58 


JR97 


36.3 


8.5 


3C 288^^ 


0.246 


31.31 


24.84 


17.32 


d96 






37.4 


8.9 


3C 293 


0.045 




23.93 


13.68 


S73 


34.84 


ZB95 


36.0 


8.9 


3C 296 


0.0237 


29.92 


23.27 


12.18 


M99 






35.5 


8.8 


3C 305 


0.0411 




23.32 


13.68 


M99 


33.41 


ZB95 


36.0 


8.7 


30 310=^ 


0.054 


30.65 


23.99 


15.66 


M99 






36.6 


8.0 


3C314.P 


0.1197 


<30.79 


<22.77 


17.91 


d96 






36.6 


7.8 


30 315=^ 


0.1083 




24.86 


15.84^= 


LL84 






36.8 


8.7 


3C317 


0.0345 


30.70 


24.29 


11.71 


GOO 


34.94 


ZB95 


36.2 


9.5 


3C 338 


0.0304 


30.61 


23.62 


12.20 


S73 


34.08 


ZB95 


36.2 


9.1 


3C 346=- 


0.162 


32.45 


25.36 


16.98 


d96 


35.35 


ZB95 


36.9 


8.8 


3L, 34o 


0. 154 


30.94 


23.98 


16.13 


E85 


33.73 


T93 


38.1 


8.9 


3C386 


0.017 


32.24 


22.24 


12.21 


M99 


<33.45 


S96 


35.4 


8.5 


3C 424=^ 


0.127 


<31.04 


24.07 


17.05'" 


ZB95 






36.7 


8.3 


3C 433=^ 


0.1016 




23.33 


14.82<^ 


L85 


35.34 


ZB95 


37.1 


9.1 


3C 438=^ 


0.29 


<31.21 


24.73 


18.24 


d96 


35.17 


ZB95 


37.8 


8.6 


3C 442=^ 


0.0263 


29.44 


21.77 


14.07 


M99 


33.85 


ZB95 


35.7 


8.0 


3C449 


0.0171 


30.36 


22.66 


13.05 


M99 


32.87 


K02 


35.2 


8.0 


3C 465 


0.0302 


30.88 


24.02 


13.14 


S73 


33.23 


M96 


36.1 


8.6 



Notes for the table 1 . a: the sources for which the jet power is greater than the maximal jet power extracted from a rapidly spinning black hole with 
g = .95 and / = 1; b: the ij-band magnitude is conveiled from the y-band magnitude assuming V — R = 0.61 (Fukusita. Shimasaku & Ichi kawal 
Il995t) : c: the i?-band magnitude is converted from the K-hwA magnitude assuming R — K = 2.5 iDunlop et alil2003 ): Column (11: source name. 
Column (2): redshift, Column (3): log^Q optical core luminosity (W A~^), Column (4): log^Q radio core luminosity at 5 GHz (W Hz-^), Column (5): 
rest frame K-band magnitude, Column (6): references for R-band magnitude. Column (7): logjg luminosity of nan'ow-lines Ha+[NlI] (W), Column 
(8): references for the line emission. Column (9): logjQ jet power estimated from total low-frequency radio luminosity (W), Column (8): the derived 
black hole masse s 

Refere nces: d96: Ide Koff et alj |l996l) . E85: Bales Igald USSSi), G77: iGodwin et alj \WT% . GOO: iGovoni et alj <200cl) . JR97: Ijackson & Rawling^ 
Il997t). LL84: citetll84 L85 : Xillv Longair & Miller (1985), M96: Mulchaev. Wilson. & Tsvetanovel 1 19961, M99: Martel, Baum & Sparks 
Martel. Baum & SDarkJll999t) . S7 3 : .Sandaga 1 1 973 ). S96 : .Simpson et aL U 996.) . T93: .Tadhunter et al-a993.) . V99: IVerdoes Kleiin et al]ll999t) . V02: 
Verdoes Kleiin et alJ<20O2ft . ZB95: lzrrbel & Baumlll995t) . 



in the sample is greater than the maximal jet power calculated for 
a = 0.95 (13 of 33 sources, hereafter referred to as high-jet-power 
sources, and the remainder as low-jet-power sources). For some FR 
I sources, their jet powers can be as high as O.OlZ/Edd- It indicates 
that if the Blandford-Znajek mechanism is the only one in opera- 
tion then it would be unable to power the jets in these sources. In 
the case of / = 10, the relation between the jet power (Jjct and 
black hole Alhh is plotted in Fig.|2| The jet power (Jjot is about 30 
times that for f — 1, which would greatly strengthen the conclu- 
sion that the Blandford-Znajek mechanism is unable to produce the 
jets in the high-jet-power sources. Furthermore, it would suggest 
that the mechanism may have difficulty powering any of the 3CR 
FR I sources. Anyway, hereafter we use, conservatively, / = 1 in 
the estimates of jet power Qjet (values for the case / 7^ 1 can be 



evaluated by multiplying by f '^^^). The uncertainty in / is there- 
fore much greater than any uncertainty in rhcrit unles s a ~ 1 
which seems inconsistent with the MHD simulations iArmitag^ 
[l998; Hawlev & Balbus 2002). 

One possible explanation may be that the jets in these radio 
galaxies are not primarily accelerated by the magnetic field lines 
threading the spinning black holes, for exampl e, they may be accel- 
erated by the field lines of the accretion disc t' Blandford & Pavna 
1982), or by the combination of the two mechanisms. It has in fact 
been proposed that the jet is accelerated more efficiently by the 
field lines threading the accretion disc (Livio, Ogilvie & Pringle 
1999; but see also more detailed calculations in Cao (2002)). Even 
if this were the case, it implies that almost all released gravitational 
energy of the accretion matter is carried away by the jets in those 
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Figure 1. The black hole mass Afbh versus jet power Qjct(assuming 
/ = 1). The solid lines represent the maximal jet power extracted from the 
spinning black holes as functions of M^jj for different labelled values of 
spin a. The sources above the uppermost solid line (for a = 0.95) are plot- 
ted as squares, noted in Table 1, and called high-jet-power sources through- 
out this paper. The filled symbols represent the sources without detected 
central compact cores (see Table 1). The symbols with a cross represent 
the sources with only upper limits on optical core fluxes. The 3CR source 
names of the high-jet-power sources are labelled in the figure, together with 
the objects 'M87' and '3C386' discussed in the text. 



sources with Qjct ~ O.OlI/Edd for f = 1. For a high value of /, 
for example, / — 10, the jet power can be as high as ~ 0.3I/Edd 
for some sources. It is hard to see how the accretion discs in these 
sources would remain in an ADAF state. However, we cannot eas- 
ily rule out the possibility that ADIOS flows are present in these 
sources, if the accretion rates of the flows are close to the Eddington 
value and the flows extend to large radii ( Blandford & Begelman 
[l999). In this case, most of the mass and energy of the accreting 
gas may be carried away in the jet, and any optical emission from 
the accretion flow could be very faint. 

Another possibility is that the accretion discs in these high- 
jet-power sources are not in ADAF state. Their accretion rates are 
in fact higher than the critical accretion rate rficrit. In this case, 
the accretion discs are expected to be the standard optically thick 
ones, and their jet power can be higher than the maximal jet power 
expected to be extracted by the Blandford-Znajek mechanism for 
ADAF cases. We note from Figs. 2 & 3 that the one object in our 
sample with a broad-line nucleus (3C386) is not one of the extreme 
high-jet-power sources, implying that standard accretion discs in 
high-jet-power FR I sources may be common in our sample, al- 
though we would need a mechanism to explain why they cannot 
easily be seen in the HST optical images. 



Figure 2. The same as Fig.Q but with / = 10 adopted in the estimate of 
Qjct- 

5.2 Optical core properties as a function of redshift 

The high-jet-power sources will also have significantly higher red- 
shifts than the other sources because jet power is estimated from 
radio luminosity and radio luminosity correlates strongly with red- 
shift in flux-limited samples like 3CR (see Fig.|3j. The median red- 
shift of the high-jet-power sources is 0.127, compared to 0.0235 
for the low-jet-power sources, so the HST linear size limits on any 
optical core emission are ~ 5 times larger for the high-jet-power 
sources (see Fig.|4}- Typically, therefore in low-jet-power sources, 
any compact optical cores have sizes < 0.033 kpc, whereas the 
corresponding limit for high-jet-power sources is < 0.15 kpc. 

Note, also, that the objects with resolved optical cores (and 
hence only limits on their optical core emission in Table 1; See 
Sec. 2) are concentrated amongst the high-jet-power sources, as is 
3C315 (see Sec. 2), the object with an elongated optical structure 
and no optical core. It seems like a large fraction of the high-jet- 
power, high-redshift sources have weak or absent cores, and evi- 
dence for optical emission extended on ~0.15 kpc scales. It is not 
clear whether the extended optical cores of these objects is associ- 
ated with an extended synchrotron component, or with host galaxy 
light. However, regardless, we can conclude that the true optical 
cores of a large-fraction of the high-jet-power sources may well be 
being obscured. 

5.3 Correlation between radio core luminosity and total 
radio luminosity 

In Fig. |5| we plot the relation between VLA radio core luminos- 
ity at 5 GHz and the total radio luminosity at 151 MHz. We find 
that a significant correlation between I/c,5G and Lisi, except two 
sources 3C28 and 3C314. 1 have relatively lower core luminos- 
ity. The generalized Kendall's t test (ASURV, Feigelson & Nel- 
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Figure 3. The total radio luminosity at 151 MHz versus redshift. 



Figure 4. The optical core luminosity versus redshift. 



son 1985) shows that the correlation is significant at the 99.98% 
level for the whole sample (the significant level becomes >99.99% 
for the sample subtracting two sources 3C28 and 3C314.1). The 
jet power estimated from the low-frequency radio luminosity rep- 
resents the power of the jet averaged over a long period of time 
(> 10^ years), while the VLA core emission (typically correspond- 
ing to sub-arcsecond angular size) reflects the emission from the jet 
near the central engine which can obviously vary at a much shorter 
timescale (~ 10''^* years). The correlation between Lc,5G and 
Li5i indicates the jet power estimated from low-frequency radio 
luminosity is a reliable measure of the current jet power in most if 
not all the cases (e.g. 3C28 and 314.1). The key point here is that 
with these possible exceptions (which could, for example, be tem- 
porarily turned-off jets) the high-jet-power sources probably have 
high jet powers throughout their life times, and the physical model 
for how they are powered must account for this. 

5.4 Optical core luminosity versus blacli hole mass 

We plot the optical core luminosity Lc,opt versus black hole mass 
Alhh in Fig.|6| The observed optical core emission may be a com- 
bination of the synchrotron emission from the jet and the emission 
from the disc, if the emission from the disc has not been obscured. 
The optical core luminosity therefore sets an upper limit on the 
emission from the disc for unobscured cases. We conve rt the optical 
core luminosity to a bolometric luminosity by using iKaspi et'all 

hoodt 



9XL 



A ,opt ■) 



(10) 



which is valid for normal bright AGNs. We plot the line c orre- 
sponding to Lboi = 0.006I/Edd iGhisellini & Celottil i200ll) pro- 
posed a diagonal line iion ~ 6 x 10~''LEdd for the division of FR 
1 and FR 11 galaxies). It is found that all sources are, as expected, 
well below this line, which indicates that most of these sources are 



accreting at very low accretion rates, assuming of course that the 
disc emission has not been obscured. 

We note from the location of 3C386 in this plot, that some 
optically unobscured objects with standard accretion discs, and 
moderately high accretion rates, exist in the FR I population. We 
have argued in Sect. 5.2 that, if we could account properly for ob- 
scuration, a significant fraction, perhaps all, of the high-jet-power 
sources might well lie in a similar region of this plot. However, the 
strong correlation between the radio and optical core luminosities 
for the low-jet-power objects, and the restrictive limits the unre- 
solved optical core places on any obscuring material (e.g. < 6pc 
for M87), together imply that at least some of the 3CR FR Is have 
genuinely low accretion rates, and are unobscured. 



5.5 Jet power versus optical core luminosity 

The relation between the jet power Qjot and optical core luminos- 
ity Lc.opt is plotted in Fig.0 Usually, the power in the jets of even 
radio-loud AGNs is only a fraction of their bolometric luminosity 
(e.g., Rawlings & Saunders 1991; Willott et al. 1999). We adopt a 
parameter ryjct to relate the jet power Qjct to its bolometric lumi- 
nosity Lboi 



Qjct — Tj jet Lhol- 



(11) 



Using Eqs. {ToJ and {TT), we can calculate the relation between 
the optical core luminosity and the jet power for a given efficiency 
factor 77jct. The lines for different values of r^jct are plotted in Fig. 

Note, first, the location of 3C386, the one object in our sample 
where we definitely have an unobscured view of the nucleus. This 
lies in a similar region of the plot to FR 11 radio sources (e.g. Willott 
et al. 1999). The location of the high-jet-power sources at high val- 
ues of jyjet 10 is explicable if, as we have postulated in Sec. 5.2, 
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Figure 5. The radio core luminosity at 5 GHz versus total radio luminosity 
at 151 MHz for the sources in the sample. 



Figure 7. The jet power, calculated assuming / = 1, versus the optical 
core luminosity. The lines con'espond to jjjct =0.01, 0.1 and 1, respectively 
where rjjct is defined by Eguation ll II FR II radio sources would lie in the 
region defined by the three parallel solid lines (e.g. Willott et al. 1999) 
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Figure 6. The black hole mass M^h versus optical core luminosity Lc,opt- 
The line represents Lbol = 0.006 L^dd- 



they have obscured nuc lei. The low-jet-power sources, which ar- 
guably are not obscured iChiaberge. Capett i, & Celotti 1999) seem 
to lie reasonably close to ?7jet ~ 1, which may be indicative of 
some difference between these objects and the population of radio 
sources with more powerful jets. 

5.6 Jet power versus line luminosity 

In Fig.|8|we plot the relation between jet power and the narrow-line 
luminosity. The generalized Kendall's r test (ASURV, Feigelson & 
Nelson 1985) shows that the correlation is significant at the 99.99% 
level for the whole sample. The linear regression by parametric EM 
Algorithm (ASURV) gives 

logio inline = 1.00(±0.16) logio Qjct - 2.02(±5.69), (12) 

for the whole sample but 3C348 is excluded (the dotted line in Fig. 
HJ, and 

logio -^linc = 1.03(±0.32) logio Qict - 3.22(±11.47), (13) 

for the low-jet-power sources (the solid line in Fig.|8}, respectively. 
It is found that the high-jet-power FR 1 sources, except 3C348, 
lie along the same correlation defined by the low-jet-power FR 1 
sources, which may imply that the narrow emission-line formation 
mechanisms are similar for both these two classes of the sources. 

5.7 Black hole mass versus line luminosity Lune 

In Fig. |9| we compare the different behaviours of the high-jet- 
power FRIs and the low-jet-power FRIs in the narrow-line lu- 
minosity versus black hole mass plane. The mean value of the 
ratio of the narrow-line luminosity to black hole mass is < 
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Figure 8. The relation between jet power and narrow-line luminosity. The Figure 9. The black hole mass versus the narrow-line luminosity plane, 

dotted line represents the linear regression for the whole sample but ex- 
cluding 3C348, while the solid hne represents the linear regression for the 
low-jet-power sources. 



logioLiino/Mbh( W Mq-I) >= 26.10 ± 0.22 for the high- 
jet-power sources, while it is < logj^Q I/iinc/-Albh( W Mq^) >= 
24.86 ± 0.16 for the low-jet-power sources. The separation of the 
two classes is clear, and is the direct result of the strong correla- 
tion between radio and narrow-line luminosities (Sec 5.6). There 
is no significant separation in terms of black hole masses: the 
mean values of the black hole masses are < logjQ A'/bh/M© > = 
8.47 ± 0.11 and < log-^^ A/bh/Mo >= 8.71 ± 0.08 for the high- 
jet-power and low-jet power sources respectively. A natural expla- 
nation is that the high emission line luminosity of these high-jet- 
power sources is mainly contributed by power re-processed from 
obscured optical nuclei. We therefore reach the same conclusion 
as Section 5.2 by independent means: high-jet-power FRI sources 
have relatively bright (i.e. as bright as 3C386), but obscured nuclei. 



5.8 Size of the torus 

The source 3C386 with an unobscured bright nucleus deviates 
significantly from the correlation between radio and optical core 
luminosities (see Fig. 5 in Chiaberee. Cacetti, & Celotti ( 1999)). 
Most high-jet-power FRI sources in our sample have higher line 
luminosities than 3C386 (see Fig. |9}, and one can infer that 
these high-jet-power FRI sources have similar or brighter intrin- 
sic optical nuclei than 3C386. The radio core luminosities of 
these high-jet-power FRI sources are in the range of ~ lO'^" — 
10"^^ W Hz^^. If the intrinsic bright nuclei were unobscured 
in these sources, as is the case for 3C386, their locations in 
Lr - Lo plane (Fig. 5 of Chiabe rge. Capetti. & C elotti ( 1999)) 
would have definitely broken down the linear correlation derived 
bv lchiaberge. Capetti. & Celotti. (119991) . This suggests that only the 



synchrotron radio and optical emission from the jets are detected in 
these high-jet-power FRI sources. 

In the receding torus model, the inner radius of 
the torus is roughly at the dust evaporation radius 
rfTevap) ~ 0.06(Lbo i /10^^ W)^''^ pc iNet zer & Laol Il993l 
iHill. Goodrich & DeFo4 Il996h . For NGC4151, it predicts an 
evaporation radius r(Tevap) — 0.015 pc (Lboi — 6.5 x 10^^ 
W, Kaspi et al., 2000), which is roughly consistent with the inner 
radius of the dust torus ~ 0.04 pc measured fro m a delay of the 
near-infrared light curve to the optical light curve iMinezaki et alj 
The evaporation radius of the dust torus in 3C386 is 
estimated to be ~ 0.02 pc (Lboi ~ 10'^^ W). The relative thickness 
of the torus s/r ~ 2 — 3 is derived from the fraction of quasars in 
the 3CR sample (Hill. Goodrich & DePov 1996). We can estimate 
the total thickness 2s of the dust tori in these high-jet-power FRI 
sources with obscured bright nuclei to be ~ 0.05 pc, because 
the intrinsic nuclei in these sources are believed to be as bright 
or brighter than 3C386. If the jet is inclined at a large angle, for 
example, ~ 45°, with respect to the line of sight, the emission 
from the accretion disc should be obscured by the torus while 
the synchrotron radio and optical emission from the jet outside 
the torus can still be detected by VLBI at pc scales or by HST at 
~ 0.01 — 0.1 kpc scales. The emission at high energy bands (e.g., 
hard X-ray bands) from the nucleus may penetrate the obscuring 
material in the torus, and then the intrinsic properties of the nucleus 
can be explored directly from its emission at high energy bands. 
However, this is applicable only if the resolution is so high (e.g., 
~ 0.05 pc) that the emission from the jet outside the torus can be 
subtracted from the emission of the nucleus. 
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6 SUMMARY 

If all 3CR FRI radio sources were powered by central engines op- 
erating in an ADAF mode, then the Blandford-Znajek mechanism 
would have insufficient power to generate the high radio luminosi- 
ties of the most powerful FRI sources. Although it is possible to 
postulate models in which an ADIOS mode is maintained in such 
sources, this would still imply that, in some cases, a significant frac- 
tion of the Eddington luminosity of a black hole can be channelled 
into jets, with negligible emission at other wavebands. 

A less contrived solution to the existence of 'high-jet-power' 
FRI sources is to postulate that, in these FRI sources at least, there 
are standard accretion discs. With the notable exception of 3C386, 
these accretion discs seem typically to be hidden (see e.g. Grimes 
et al. 2003 for discussion of why this might be the case). The large 
fraction of the high-jet-power FRI sources that lack any unresolved 
HST optical core is consistent with this idea, as in these cases the 
true nucleus is presumably hidden, and any optical synchrotron is 
coming from scales larger than the obscuring material. Also con- 
sistent with this idea is the strong correlation between narrow-line 
and radio luminosities which means that, at fixed black hole mass, 
high-jet-power FRI sources have intrinsically brighter optical cores 
than low-jet-power FRI sources, but these cores must be hidden 

If the efficiency factor / (see Equation |2j is nearer to ~ 10 
than ~ 1 (as suggested by Blundell & Rawlings 2000) then these 
same arguments can probably be extended to include the majority 
of 3CR FRI radio galaxies. Perhaps then, essentially all of the 3CR 
FRIs have standard accretion discs. The case of M87 provides the 
best example of an object with a very low accretion rate, and no 
torus bigger than ~ 6 pc. However, these two factors alone are not 
sufficient to prove that M87 is accreting in a regime described by 
an ADAF or ADIOS model. Indeed the location of M87 in Fig. 3, 
suggests that, if / ~ 10, then its jets may also be too powerful 
for it to be plausibly associated with an ADAF mode of accretion 
unless there are powering mechanisms which dominate over the 
Blandford-Znajeck process. 

We conclude that the idea that the physical cause of the di- 
vision between FRI and FR II sources is that they have different 
modes of accretion (e.g. Ghisellini & Celotti 2001) is not supported 
by the radio and optical properties of 3CR FRI radio galaxies. This 
is in line with arguments made on the basis of individual objects 
such as the 3CR FRI 3C386 ( Simpson et al. 1996) and the 'radio- 
quiet' FRI quasar E1821+643 iBlundell & RawUngj200"ll) . 
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NOTE ADDED IN PROOF 

After the final version of this paper was submitted, we heard that 
the nucleus of 3C386 is superimposed by a foreground star (D. 
Marchesini, private communication). However, a broad Ha line 
is obviously seen in the spectrum of this source (Simpson et al., 
1996), The equivalent width of this broad line corresponding to the 
central compact core emission measured by Chiaberge et al. (1999) 
is ~ 15 A, which implies that the optical continuum emission from 
the nuclei of 3C386 is not obscured and still quite bright, though it 
may not be so bright as observed due to the contamination by the 
superimposed star. The broad line and continuum emission could 
be variable as is common in low radio luminosity broad-line radio 
galaxies (BLRGs) like 3C390.3. This may be verified by further 
optical spectroscopic observations on 3C386. 
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